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Abstract. Models for crowd dynamics are presented and compared. Well 
posedncss results allow to exhibit the existence of optimal controls in various 
situations. A new approach not based on partial differential equations is also 
briefly considered. 

1. Introduction. From a macroscopic viewpoint, a moving crowd can be described 
through its density p = p(t, x), a function of time t € M + and space i£l 2 attaing 
values in [0, 1]. In standard situations, the number of pedestrians is conserved, so 
that J R2 p(t, x) dx is independent of t. Hence, it is natural to use the conservation 
law 

d t p + div x (pV)=0. (1) 
Any model of this kind depends on the speed law that defines the velocity V of the 
crowd as a function of t, x, p, . . . A simple version of (1) is obtained assigning 



V = v(p) v(x) with 




In this case, Kruzkov Theorem [24, Theorem 1] applies and ensures that the Cauchy 
problem for (l)-(2) has a unique solution in C ' 1 (R+; L 1 (M 2 ; [0, 1])) which depends 
Lipschitz continuously from the data and, by [12, Theorem 2.6], also from v and v. 

According to (2), at time t the pedestrian at x moves along a prescribed tra- 
jectory, an integral curve of v, with a speed v(p) that depends on p evaluated at 
point x and time t. On the contrary, Section 2 is devoted to (1) with the speed 
of the individual at x depending on an average of the density p in a neighborhood 
of x. The resulting model has a rich analytical structure, the solutions being also 
differentiable with respect to the data and to the speed law. 

In Section 3 the direction chosen by the pedestrian at x depends from an average 
of the density gradient Vp around x, while his/her speed depends from p evaluated 
at x. The resulting solutions display qualitative properties usually seen in context 
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where individuals have a proper volume such as the Braess paradox [3] and the 
formation of queues [23]. 

If the various individuals have different destinations then it is possible to subdi- 
vide the crowd under consideration into different, say n, populations with densities 
pi, . . . , p n , each having a different destination. The resulting model 

d t pi + div x (pi Vi) = i = l,...,n (3) 

consists of a system of conservation laws that, when n = 1, reduces to (1). The 
results in both Section 2 and Section 3 can be extended to this more general setting. 

Finally, Section 4 approaches the problem of driving a crowd with a few moving 
individuals. First, a model based on (1) is recalled and then an approach based on 
differential inclusions is presented. The latter approach, developed following [5, 6], 
neglects the crowd internal dynamics and allows for a simpler analytical framework. 

We refer for instance to [2] for an account of the fast development of the recent 
macroscopic modeling of crowd dynamics. Moreover, measure valued conservation 
laws were considered in [18, 26]; the results in [25] deal with constrained velocity 
models; various ID attempts are found in [I, 15, 16, 20, 21]. Throughout, for the 
basic results in the theory of conservation laws we refer to [4, 19] . 

2. NonLocal Speed Choice. Consider (1) with the nonlocal speed law 

V(p) =v(p(t)*r]) v. (4) 

Here, the speed v at time t of the pedestrian at x depends on the averaged density 
(p(t) * ?y) (x) = L 2 p(t, x — y) rj(y) dy. The direction of the velocity is given by the 
(fixed) vector v{x). 

For simplicity, we state the results below in R 2 . However, the case where the 
region available to the crowd is constrained by, say, walls or doors can be easily 
recovered in the present framework, along the technique used in [7, 8] 

As is typical whenever Kruzkov techniques apply, space dimension 2 plays no 
role and the results below can be extend to R™. 

Existence and uniqueness of a solution to the Cauchy problem for (l)-(4) follow 
from the next result. 

Theorem 2.1. [9, Proposition 4.1], [10, Theorem 2.2] Let v G (C 2 nW 2 '°°)(R;R), 
v G (C 2 n W 2 ' 1 )^ 2 ^ 2 ), T) G (C 2 n W 2 '°°)(R 2 ;R). Assume Po G (L 1 n L°° n 
BV)(R 2 ;R + ). Then, (l)-(4) with initial condition p Q admits a unique weak entropy 
solution p G C° (R + ; L : (R 2 ; R + )) . Furthermore, we have the estimate LOO ^ 

IIPoIIloo 6 ^*; where the constant C depends on v, v and r\. 

The definition of weak entropy solutions is based on Kruzkov notion [24, Defini- 
tion 1], see also [9, 10]. The proof relies on a contraction argument based on the 
key estimates provided by [12, Theorem 2.6]. 

Another contraction argument, based on tools from optimal transport theory, 
allows to extend the above result to the measure valued setting in [17]. (Below, 
.M + (R + ) is the set of positive Radon measures on R 2 ). 

Theorem 2.2. [17, Theorem 1.1] Assume v G (L°° n Lip)(R;R) ; v G (L°° n 
Lip)(R 2 ;R 2 ), r\ G (L°° n Lip)(M 2 ;R+). Let Po G M + (M. 2 ). Then, there exists a 
unique weak measure valued solution p G L°°(R + ; A4 + (M. 2 )) to (l)-(4) with initial 
condition Po . Lf furthermore p Q G L 1 (R 2 ;R+), then p G C° (R+: L 1 (R 2 ; R+)) . 
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In general, in (l)-(4) no a priori uniform L°° bound on the density is possible. 
Indeed, assume that the density is 1 all along the trajectory of the pedestrian at 
x. The averaged density around x may well be less than 1, forcing the pedestrian 
to proceed and, hence, leading to a increase in the density. This behavior can be 
related to the rise of panic, see [15, 16]. In the literature, values of p of up to 10 
individuals per square meter were measured, see for instance [22]. 

Aiming at preventing the insurgence of these phenomena, it is natural to consider 
control problems where functionals of the density of the type 

Jt(Po) = / / f (p(t,x)) dx dt where p solves (l)-(4) with datum p (5) 
Jo Jn 

have to be minimized. Here, Q is the region where the density needs to be controlled 
and / is a C 1 function weighing on acceptable densities and quickly increasing 
when p approaches dangerous values. Necessary conditions for the minima of (5) 
are available once the differentiability of the solution to (l)-(4) with respect to the 
initial datum is proved. This motivates the following result. 

Theorem 2.3. [9, Theorem 4.2] [10, Theorem 2.2] Let Po G (W 2 '°°nW 2 ' 1 )(R 2 ;t+), 
r Q £ (W^nL 00 )^ 2 ;^). Assume ve (C 4 nW 2 '°°)(R; R), v G (C 3 nW 2,1 )(R 2 ; R 2 ), 
r\ G (C 3 H W 2 '°°)(R 2 ; R + ). Then, there exists a unique weak entropy solution 
r G C°(R+;L 1 (R 2 ;R)) to the Cauchy problem 

dtr + div (rv(p * rf) v{x)) — — div (pv'(p * rj) v{x)) , r(0) = r Q . (6) 

Furthermore, for all p a G (W 2 ' 1 n W 2 -°°)(R 2 ; R+) and r D G (W 1 ' 1 n L°°)(R 2 ;R), 
call ph the solution to (l)-(4) with initial datum p + hr . Then, for all t G R + , 

Ph(t)-p(t) 



lim 

h-K) 



-r(t) 



= (7) 

L 1 



i.e., the solution p to (l)-(4) is Gateaux differ entiable in p a along any direction r Q . 

To prove this theorem, first the well posedness of (6) is obtained and then the 
limit (7) is computed. In both steps, the estimates in [12] play a key role. At present, 
no analog to Theorem 2.3 is available in the setting of Theorem 2.2. Indeed, a good 
definition of Gateaux differentiability on the set of probability measures equipped 
with the Wasserstein distance of order 1 is, to our knowledge, not available. 

3. NonLocal Route Choice. Consider (1) with the nonlocal speed law 

V(p)=v(p) (v{x)+l{p)) . (8) 

Here, the individual in x at time t moves at the speed v (p(t,x)) that depends on 
the density p(t, x) evaluated at the same time t and x. The vector V{x) € R 2 is the 
preferred direction of the pedestrian at x, while X[p)(x) describes how the pedestrian 
at x deviates from the preferred direction, given that the crowd distribution is p. 
Thus, the individual at time t in x is assumed to move in the direction of the vector 
i>{x) + [x (p{t))^j (x). The basic well posedness result for (1)— (8) is the following. 
Theorem 3.1. [8, Theorem 2.1, Theorem 2.2] Let the following conditions hold: 

(v): v G C 2 (M; R) is non increasing, v(0) ~ V and v(R) = for fixed V, R > 0. 

(z7) : y g (C 2 n W 1 - 00 )(R 2 ;IR 2 ) is such that div z? G (W 1 - 1 n W 1,00 )(IR 2 ; M). 

(I): 1 G C° (I/QR 2 ; [0, R]); C 2 (M 2 ; R 2 )) satisfies the estimates: 

(1.1) There exists an increasing Ci G L^ C (R + ;R + ) such that, for all r G 
I/(R 2 ; [0,R]), ||l(r)|| I>0O < C7(||r|| L1 ) and ||divJ(r)|| Ll < C 7 (||r|| L1 ). 
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(1.2) There exists an increasing Ci € L^ C (R + ; R + ) such that, for all r G 
I/(R 2 ; [0,R]), ||VdivI(r)|| L1 < Cj(\\r\\ L1 ). 

(1.3) There exists a constant Kj such that for all r\,r 2 G L 1 (M 2 ; [0,R]), 

\\T{ ri ) -X{r2)\\^ < A-j-Hn-rallLi, 
||X(n) -X(ra)|| tl + ||div(X(ri) -X(rs»))|| Ll < Ki ■ ||n - r 2 || L1 . 

Choose any p G (L 1 fl BV)(R 2 ; [0, R]). Then, there exists a unique weak entropy 
solution p G C° (R+;L 1 (M: 2 ; [0,R])) to (l)-(8). Moreover, p satisfies the bounds 

||pW||ti = \\Po\\ L i , for a.e. t eR + , 



TV(p(t)) < TV( Po )e' £t + I ie fct 7V|| g |! LOO([0ifl]) (||Vdiv l ;|| Ll +C / (||p o || Ll )) , 



where k = (2N + l)\\q'\ 



L°°([0,.R]) 



(||Vz?|| Loo + C/(||p || L i)) . If also the speed law 



V'(p)=v'(p) (v>{x)+T'{p) 



(9) 



satisfies the same assumptions, then the solution p to (l)-(8) and p' to (l)-(9), 
with data p ,p' G (L 1 n BV)(R 2 ; [0,R]), satisfy 

\\pi{t) - P2 {t)\\ L1 < (l + Cr(t))||poA-^,a|| L i + C(t)(||? 1 -ga|| w i, 00 + <i(Ii,^ l )) 
+C(t) (Vi - z7 2 || LDO + ||div fa - z7 2 )|| L1 ) 

where 

d(T x , X 2 ) = sup 1 1 |J X (p) - 1 2 (p) 1 1 ^ + 1 1 div (1, (p) - 1 2 (p)) 1 1 ^ : p G L 1 (R 2 ; [0, R] ) } . 

The map C G C°(R + ;R + ) vanishes at t = and depends on TV(po,i) ; ||A),i|| L i> 
II^iIIloo, ||divz7i|| wlil , |jgi|| w i.oc, ||g2|| w i>°°- 

In operation research, Braess paradox states that adding extra capacity to a net- 
work can, in some cases, reduce the overall performance of the network, see [3]. A 
relevant problem in the design of escape routes is the planning of suitable devices 
that reduce the exit time. The model (l)-(8) allows to show that the careful intro- 
duction of suitable obstacles in suitable locations does indeed reduce the exit time. 
In fact, these obstacles reduce congested areas at the sides of the door jambs. 



Time= 0.000 Inside = 100.00% 




I 



v{p) 



6(1 -p), 



fa) = [l"(^)J 
Po (x) = 0.75x [27]x[ _ 22] (x), 



r = 0.6, e = 0.4. 
Figure 1. Initial datum and room geometry, without obstacles. 



(10) 



We consider a room with an exit, as in Figure 1. The vector v = i>{x) is the 
unit vector tangent at x to the geodesic connecting x to the exit and T{p) = 

1~2 



-e {V(p*rj)) I \Jl+\\V(p*v)\ 



see (10). 
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Figure 2. Solution to (l)-(8)-(10) with e = 0.2, at times t = 
4.438, 6.253, 11.396. On the first line, no obstacle is present. On 
the second line, 4 columns direct the crowd flow. The exit time in 
the latter case is shorter than in the former one, see [7]. 



The careful positioning of obstacles as in the second line of Figure 2 diminishes 
the size of the congested region and, with the chosen initial datum, gives an exit 
time lower than that with the room free from any obstacle, see Figure 2. 



4. Individuals Driving a Population. We finally introduce a model describing 
the situation in which a discrete set of isolated individuals interacts with a contin- 
uum crowd. Examples can be a (group of) predator(s) running after their preys, 
shepherd dogs driving a herd of sheep, or a leader attracting a group of followers to 
a given region. Let p £ M + be the population density and p = (pi,... ,pk) £ R 2fc 
be the positions of the k individuals. Following [11], the interaction is described by 

dtp + div (^pV(t,x, p(t,x),p(t))^ = 0, , n > 
p = tp(t,p(t),p(t)) . 

Here, ip is typically nonlocal, meaning that the individuals p react to averages of 
quantities depending on p. The well poscdness of (11) is proved in [11, Theorem 2.2], 
by means of Kruzkov theory, the estimates in [12] and tools from the stability of 
ordinary differential equations. 

As a first illustrating example, assume that the vector p £ R 2 is the position of 
a leader (e.g. a magic piper) and p is the density of the followers (e.g. rats). We are 
thus lead to consider (11) with 

V(t,x,p,p) = v(p) (p - x) e-IIP-*!! 

<p(t,*, P ) = (i+( P * v )( P (t))) m. (12) 

The function v essentially describes the speed of the followers and is, as usual, a 
smooth decreasing function vanishing at, say, p = 1. The follower located at x 
moves along p(t) — x toward the leader, with a speed exponentially decreasing with 
the distance \\p — x\\ between leader and follower. The speed of the leader increases 
with the averaged density p * r], computed at the leader's position. Indeed, we 
expect the leader to wait for the followers to join him when the followers' density 
around him is small. The direction ij) of the leader is chosen a priori. See Figure 3 
for a numerical integration of (11)— (12) and [11] for further details. 

As a further example, consider n shepherd dogs, located in Pi(t) £ IR 2 for i £ 
{1, . . . , n } and a group of sheep of density p. The dogs have to confine the sheep 



T? AA nr\j ni\;inn \a had aa/t?t to i\a t t-t , /~ i t tt? "ptt t y t\ t t? f? z -1 t t? f? w pnnnnAru 






FIGURE 3. Solution of (11)-(12), from [11]. 



within a given area. We are thus lead to consider (11) with 

V(t,x,p,p) = v{ P )v(x)+Yiu(x-p l y-^- x * ,] 

<p(t,v,p) = , (p * v "^ (p ' (t)) 2 , farie {l,...n}. 



l+|lp*Vr;(p, (*)) | 



(13) 



As above, the speed of the sheep is given by the decreasing function v(p) that 
vanishes in p = 1. The direction of a sheep located at x is a sum of two terms. 
The first one is the sheep's preferred direction u(x); the second one is the vector 



£i=l(z-Pi(*))« 



' representing the repulsive effects of the dogs on the sheep. 



Each dog runs around the flock along the direction perpendicular to the gradient 
of the sheep average density. 



5. A Different Approach. Following [6, 13, 14], we present another framework to 
describe the population-individuals interactions. Initially, the population occupies 
the compact set K Q C R 2 . If there are no individuals, the member at x of the 
population is free to wander in R 2 , according to the differential inclusion 

ieB(0,c), x{0)eK o , (14) 

c being the maximal wandering speed and B(0, c) the closed ball in R 2 centered at 
with radius c. Hence, the population fills the reachable set of (14). Introduce now 
?? individuals sited at £ = (£i, £2, ■ • ■ , £n) G R 2 ™. Then, the interaction between the 
individuals and each population member leads to the modified differential inclusion 

x€v(x,£(t))+B{0,c), x(0)&K o , (15) 

where the vector field v € C 0,1 (R 2 x R 2 ™; R 2 ) is the drift speed due to the attractive 
or repulsive effect that each agent has on each member of the population. Thus, 
given the individuals' trajectory £ € C°^(R + ;R 2n ), the reachable set TZ^(K ,t) 
of (15) at time t is the set occupied by the population at time t under the effect of 
the agents. With the present assumptions, Tt^(K a ,t) is non-empty and compact. 
If only one agent is present (n = 1) and v is spherically symmetric, i.e., 

v(x,£) = ip(\x - £\)(x - for a suitable tp: K -> R. (16) 



the next result exhibits a trajectory £ confining the population in a given set K. 
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Theorem 5.1. [14, Theorem 2.8] Let c> 0. Fix a bounded ip G C^(R n ;R) and 
define v as in (16). Assume that there exist positive R~ , Rj and R such that 

~ jf V [y R 2 + Rl - 2R* R cos6^ (R* - R cosfl) d(9 < -c 

/or a// i?» G [iZ~ , . Then, there exists o(€ C° * (R+;dB{0,R)) such that, 
calling K a the region initially occupied by the population, 

if K C B(0,R~) then H^(t,K ) C B(0,R+) for all t>0. 

Note that the confining strategy t n- £(t) above is constructed explicitly, see [13, 
Theorem 2.5]. A negative result is also available. Before stating it, recall that 
for a measurable function tp : R + — > R, its non- decreasing rearrangement is the 
function <p*: R + — > R, which is non-decreasing and satisfies C 1 (p* 1 (]— oo, a])) = 
/I 1 (yj-^J-oo,^)) for all aeR. 

Theorem 5.2. [14, Theorem 2.7] Lef c > 0. Fix a bounded ip G C^(R;R) and 
define v as in (16). Le£ (y9* 6e the non-decreasing rearrangement of the function 

If the initial set K Q is such that 

2cy/^b z + f V3*(s)ds>0 for all a>C 2 {K ) 
Jo 

then, for every £ G C^(R + ;R fc ™) ; the measure C 2 (TZ^(t, K a )) of the reachable set 
TZ^(t, K a ) of (15) increases unboundedly in time, so that no confinement is possible. 

We refer to the cited references for the statement of these results in arbitrary 
space dimension. Theorem 5.2 holds also in the case of several individuals, each 
acting as in (16) (see [14]). 
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Section at time Q.194 
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Solution at time Q.3Q3 



+ 



1 ^5 



03 



Solution at time Q.364 




1 ^5 



03 



Solution at time 0.400 
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